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Formation of localized magnetic states in a metallic host is a classic problem of condensed matter physics 
formalized by P. W. Anderson within the so called single impurity Anderson model (SIAM). The general picture 
in a host of a simple one-band metal is that a large Hubbard U in the impurity orbital is pre-requisite for the 
formation of localized magnetic states. In recent years three dimensional (3D) Dirac solids have emerged the 
hallmark of which is strong spin-orbit interaction. In this work we show that such a strong spin-orbit interaction 
allows to form localized magnetic states even with small values of Hubbard U. This opens up the fascinating 
possibility of forming magnetic states with s or p orbital impurities - a different from traditional paradigms of 
dot f orbital based magnetic moments. 
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I. INTRODUCTION 


Dirac matrices, namely: 


When an element possessing usually dov f orbital is added 
as an impurity to a metal, under certain conditions the impu¬ 
rity atom can have a magnetic moment. This problem was 
formulated and solved within a mean field approximation by 
AndersorP. The essential ingredient was identified to be the 
Hubbard U which is relatively large in d and / orbital impuri¬ 
ties. Anderson’s formulation lead to the following picture: If 
the Hubbard U is large enough double-occupancy and empty 
charge configurations of the impurity become energetically 
costly. If the hybridization V with continuum of Bloch states 
in the host metal is weak enough to prevent decay of the local¬ 
ized spin-split states into continuum, within the Haitree mean 
field it leads to the formation of localized magnetic states in 
metallic hosts. In a simple metal considered in the original 
Anderson impurity model the spin-orbit is absent and there¬ 
fore the study of interplay between the spin-orbit interaction 
and other parameters of the Anderson impurity problem re¬ 
mains an outstanding problem. 

One interesting paradigm where strong spin-orbit interac¬ 
tion manifests itself in a fascinating way is concerned with 
three dimensional Dirac solids. Dirac electrons in solids 
appear under quite general conditions where in presence of 
strong spin-orbit interactions, two bands of Kramers doublets 
are separated by a small gajM Under such general conditions 
the effective bands of the solid can be represented by the Wolff 
Hamiltonian: 


iT-vv = A 7 ° -F k. ^ 
i=i 


( 1 ) 


where v is related to veloci^ matrix elements and 7 ^ and 7 ° 
are 4 X 4 matrices given b 


:iw 1 


7 ° = 


I 0 
0 -I 


7-^=1 


0 
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( 2 ) 


with I and cr^ ,j = 1, 2, 3 being unit and Pauli 2x2 matrices. 
The four 7 ^ matrices with ^ = 0,1, 2, 3 satisfy the algebra of 


= 2 <'^ 


(3) 


where the matrix representation of the tensor in the right side 
is 77 ^“^ = diag(l, —1, —1, —1). Note that this representation 
of Dirac matrices is slightly different from those used in the 
field theory textJ^. 

The isotropic Wolff Hamiltonian corresponds to situation 
where velocity matrices are the same for three Cartesian di¬ 
rections: Vj = vCj with Cj being three mutually orthogonal 
unit vectors along the x, y, z directions. Under isotropic con¬ 
ditions the Wolff Hamiltonian becomes. 


Hjy = 


A 

—ivk.cr 


ivk.a \ 

-A J- 


(4) 


This is precisely the Dirac Hamiltonian where the light ve¬ 
locity is replaced by v. The above effective Hamiltonian is 
obtained by k.p expansion around a particular wave-vector kg 
(corresponding to L point in bismuth). 

One of the long standing puzzles of such Dirac systems 
in the context of bismuth has been the Diamagnetic response 
that was markedly different from normal metallic states. Un¬ 
like the one-band situation where a Landau-Peierls formula 
describes the diamagnetic response of the solid at hand, in 
the case of Dirac systems the inter-band effects play a crucial 
rolJ^. In this context an interesting question can be formu¬ 
lated with respect to the behavior of impurity states in a host 
of Dirac electrons: What is the role played by the presence of 
the other band of Kramers doublets? Another interesting as¬ 
pect of local moment formation in 3D Dirac solids is related 
to the role of spin-orbit interactions that determines the veloc¬ 
ity scale V in this Hamiltonian. What is the interplay between 
the spin-orbit interaction encoded in v and the formation of 
local magnetic states? In this work we will show that, un¬ 
like normal metals where basically a strong Hubbard U at the 
impurity orbital causes local magnetic moment formation by 
excluding double occupancy, in the case of 3D Dirac solids 
the spin-orbit interaction facilitates the formation of localized 
magnetic states even with very small values of Hubbard pa¬ 
rameter U. This has far reaching consequences: In addition 
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to impurity atoms with d or f orbital, even systems with s or 
p impurity orbitals having smaller values of Hubbard U have 
a chance of forming magnetic moments in 3D Dirac solids - 
an opportunity not available in a host of normal metal. More¬ 
over the presence of two bands in a 3D Dirac solid produces 
an additional region in the phase diagram which is which has 
no counterpart when the host is a normal metal or even a 2D 
Dirac system (graphene). We clarify that this new portion of 
phase diagram can be considered as the signature of a second 
band of Kramers doublets. 

The paper is organized as follows: In section II we formu¬ 
late the single impurity Anderson model (SIAM) in a host of 
3D Dirac solid. We set up mean filed equations parameteriz¬ 
ing hybridization of impurity orbital with local orbitals of the 
2D Dirac host from a purely A-sublattice type to B-sublattice 
type. In section III we report our numerical results leading 
to section IV on applications to realistic Dirac materials, in 
Section V we discuss possible deviations from a simple Dirac 
Hamiltonian such as the presence of anisotropy as in the orig¬ 
inal Wolff Hamiltonian or addition of a quadratic term that 
could possibly give rise to a non-trivial topology of host. 


II. FORMULATION OF THE PROBLEM 


The isotropic Wolff Hamiltonian for a general 3D Dirac ma¬ 
terial is given bjEl 




A ivk.d 
—ivk.S —A 




fc- 


(5) 


To construct the hybridization part, consider the elemen¬ 
tal bismuth corresponding to the Hamiltonian (|^. Bi has 
a rhombohedral lattice with bases composed of two atoms. 
The creation operator d) ) creates states in symmetric (anti¬ 
symmetric) orbitals corresponding to the top of valence band 
(bottom of conduction band). Therefore if we define a^{b^) as 
creation operator on sub-lattice A{B) then ~ -f and 
fd ^ (jt _ Therefore a local hybridization with an orbital 
on site A or i? is given as ± Here we assume a 

quite general combination of c and d states as Ac -I- ^d. There¬ 
fore the hybridization of impurity with the host electrons is 
assumed to be: 

k 

(9) 

where V is hybridization strength between impurity level and 
the Bloch states. As an example one can set A = 1), ^ = 

-fl to hybridize the impurity with just one atom from sub¬ 
lattice A{B). On the other hand A = -1-1, ^ = 0 represents the 
anti-symmetric hybridization with two sub-lattices (conduc¬ 
tion), while A = 0,^ = -fl stands for symmetric hybridiza¬ 
tion with two sub-lattices (valence). If one writes the equation 
of motion for the impurity’s Green function in frequency do¬ 
main, aUlfl))) it gives: 


(w — £o){{fa\fl')) — + 

U{{fan-Ml)) + ^ ^((e,-J/],)), (10) 

k 


with a four component spinor T't = ( ct , ct , dt , dt ). 

In the above basis the operator c^^(dj^^) is creation opera¬ 
tor in the upper (lower) band. These operators at k — ko 
(corresponding to L point in bismuth) can be thought of as 
anti-bonding (bonding) molecular orbitals composed of two 
atomic orbitals at two locations of a unit cell. This Hamilto¬ 
nian represents a gapped two bands model where 2A is the 
energy gap and v is the velocity of Dirac fermions which is 
usually 2-3 orders of magnitude smaller than the light veloc¬ 
ity. The spectrum of this Hamiltonian is : 

Efc = ±\/ + A^. ( 6 ) 

As the canonical model for the study of magnetic moment 
formation in the host of itinerant electrons, we consider the 
SIAM as follows: 


where we define = Acg^ + ^d^^. The Hartree approx¬ 
imation in this case corresponds to replacement Eq ^ = 

^0 + U{ng). Writing the equation of motion for ((eg^l/j^/)) 
yields: 

{((M/^^))[(w-A)|Ap + (a; + A)|e|2+*a(Ar-A*a^'fc.] 

+ - laky)]. (11) 

By combining Eqns. ( [T0| and ( [TT] i, the impurity’s Green func¬ 
tion is given by: 

(w — Co- — Ej(tLi))((/£r|/j^;)) = Saa', (12) 


H — Uq ddinip ddjiyb, 


(7) 


where the self-energy is given by: 


where impurity’s contribution is: 

TTimp — E e-dflfcr + Un^n^. (8) 

(7 

Here represents creation operator in impurity level of en¬ 
ergy £o- U is the Hubbard repulsive interaction in impu¬ 
rity site. 




|U|2 (|A|2 + |g|2) + A(|^|2-|A|2) 

Af ^ ixP- — el 


In order to obtain the above self-energy the fact that if a 
given state at k is occupied, the time-reversed state at —A: is 
occupied too simplifies the integration. Summation over k 
for obtaining self-energy leads to diagonality with respect to 
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Hybridization of impurity with conduction band (A, = 1, ^=0) 



X = D/t)^ 
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site, i.e. and n^. Occupation number can be calculated as 
follows; 


^)=-- 


duj- 




[Z-i(w)w-es] +5(S/)2 


(15) 


where. 


2ttujD 


' , , y/D^ - A2 + - A2 ^ Vi92 - A2 

21n(^^ - / ) — 4 , 

\ /-■ A- ,.■ A- Va;2 - A2 


(16) 


Difference of occupation number is local magnetic moment, 

m = — n^. 


III. NUMERICAL RESULTS 


FIG. 1. (Color online) Phase diagram for local moment formation 
in 3D Dirac solids. The area enclosed by the curve and Y -axis is 
magnetic region. The subscript fj, in V), emphasizes that for obtaining 
this diagram chemical potential has been tuned. In this case, impurity 
level lies in conduction band. 


spin and simplifications in the Green function for general hy¬ 
bridization pattern parameterized by arbitrary A and Per¬ 
forming the integration over k in Eq. ([T3), the self-energy as 
a function of uj becomes; 


, ^(iAp + m+A(i;p-iAn .^^ 

2ttujD 


VD2 - A2 -F Vw2 - A2 - A2 . 

^ y'D'^ - A2 - Vw2 - A2 ^ A/a;2 - A2 ' ' 


(14) 


In the above equation D is the bandwidth cut-off and the 
parameter v = V'^D/{2 ttv^) naturally emerges in the the¬ 
ory that in addition to the Fermi golden rule decay rate pro¬ 
portional to hybridization strength V contains the velocity 
scale V (that arises from the spin-orbit coupling). This means 
that in a host of 3D Dirac electron the hybridization strength 
and spin-orbit coupling do not independently determine the 
physics of local moment formation; instead the specific com¬ 
bination ~ V'^ I plays the role played by the combination 
V'^/U in normal metals. This is already very suggestive that 
in 3D Dirac solids a large spin-orbit coupling may lead to 
local magnetic moments in a similar way the Hubbard U in 
normal metals does. As will be seen in the next section, this is 
indeed the case, and unlike normal metals, the 3D Dirac solids 
allow for magnetic moment formation even for impurity or¬ 
bitals with small values of Hubbard U. We therefore use the 
combination X = V^Djv‘^ to construct the phase diagram of 
local magnetic moment formation in 3D Dirac materials. In 
order to contrast the phase diagram against normal metallic 
hosts, we choose Y — {p, — eo)/U. Moreover, Eq. ( [I4| ) shows 
that the relative phase of complex numbers A and is not im¬ 
portant in phase diagram which is determined by self-energy. 
We compute occupation number for both spins in impurity’s 


We solve Eq. for t and I spins self-consistently. By 
tuning V (X axis) and p (which results in scanning Y axis) 
we construct the region of parameter space corresponding to 
localized magnetic states in three dimensional Dirac materi¬ 
als. We consider the case of impurity level with eg = 40 meV 
and U = 1 eV which hybridizes anti-symmetrically with two 
sub-lattices; this leads to hybridization with conduction band. 
Phase diagram is presented in Fig{^ This result is suggesting 
that stronger spin-orbit coupling which means larger velocity 
V of Dirac fermions, leads to larger magnetic moments (red 
color in the intensity plot means magnetic moments closer to 
one - in units of h/2). Due to spin-orbit coupling, the true 
eigen-states of the Hamiltonian Hg of the host Dirac material 
is actually a Kramers doublet which is a linear combination 
of spin t and | states. For large enough U (small Y) where 
conditions for single-occupancy of the impurity orbital are fa¬ 
vorable, the mechanism that can reduce the local moment e.g. 
in the t state is the tunneling out of the impurity f state to an¬ 
other t state in the surrounding continuum of states of the host 
material. But since the eigen-states of the Hg are not purely 
t nor purely states, the spin-orbit coupling weakens the rate 
of such transitions out of and into the localized state. There¬ 
fore stronger spin-orbit coupling is expected to give rise to 
stronger localized magnetic moments under comparable con¬ 
ditions. The formal way of understanding the above situation 
is that the broadening of spin-split impurity states are given 
by r ^ V^p{e)/v^. Therefore larger spin-orbit coupling (v) 
leads to smaller broadening, and hence a more perfect spin- 
split levels, i.e. larger magnetic moments. This form of spin- 
orbit coupling appears only in 3D Dirac solids, and hence the 
present mechanism of the local moment enhancement can be 
considered as a characteristic of these systems. 

The second property of the region of magnetic moments in 
the above hgure is that it is conhned to F >0. While this fea¬ 
ture is similar to the behavior of localized magnetic moments 
in a host of normal metallic hosts, it is distinct from the mag¬ 
netic moment formation in the two dimensional Dirac systems 
(graphenej^. As pointed out in Ref.0 this difference can be 
traced back to the damping behavior of the local Green’s func¬ 
tion at large ui. This behaves as in normal metals and 3D 
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Hybridization of impurity with conduction band (A, = 1, ^=0) 
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FIG. 2. (Color online) Phase diagram for local moment formation 
in 3D Dirac solids. The impurity level £o = —1.4 eV lies deep in 
the valence band. The color code is indicated in the legend. The 
subscript /r in means that to obtain this diagram chemical poten¬ 
tial has been tuned. In this figure the tail of the magnetic phase is 
extended along ^ — eq ~ U line. 


Dirac solids, while in the 2D Dirac systems it goes as 
Despite that local magnetic moment formation in 3D Dirac 
and normal metals both occur for ^ > 0 (i.e. > Eq). it 

is interesting to note that the upper boundary of the magnetic 
region has positive curvature in 3D Dirac solids, while in the 
normal metals the curvature is negative and the upper bound 
is convex. 

In Fig. 1^ we construct the phase diagram for a very nega¬ 
tive value of £o = —1.4 eV. This figures shares the general 
property with Fig. [^that larger spin-orbit coupling gives rise 
to stronger magnetic moments. However they differ in their 
large X behavior. In Fig. [^corresponding to Sq = 40 meV, 
weak magnetic moments for large X values are formed when 
/i — Eo « 0+, while in Fig. [^ corresponding to Eq = —1.4 
eV, the corresponding small moment states are formed for 

— Eq ^ U. The large X is equivalent to small spin-orbit 
coupling and large hybridization V. In the limit of large X, 
the spin-orbit coupling becomes negligible. Assuming that 
the Hubbard U is large enough to favor single occupancy, the 
dominant term to be minimized will be hybridization term giv¬ 
ing rise to energy contribution \V\'^p{^). For a given Eq and 
allowing p, to be variable, the minimization of the above en¬ 
ergy contribution amounts to selecting regions with smaller 
DOS. For small values of Eq = 40 meV in the conduction 
band subject to the y > 0 condition, the smallest value of 
p{p) is obtained when p —> Eq. In the case of very nega¬ 
tive Eq = —1.4 eV, minimization of p{p) is achieved when 
p ^ Eq + U. The sign of Eq affects the elongation pattern of 
magnetic region in the same way as graphene^. 

So far we have been concerned with the situation where 
the impurity orbital was coupled to the conduction band, i.e. 
A = +1,^ = 0. It is interesting to see what happens when 


Phase diagram for different hybridization pattern 



FIG. 3. (Color online) Comparison of different hybridization in 
which we normalize strengths. In this plot the chemical potential 
is tuned and the impurity level and Hubbard U set eq = 40meV and 
U — leV respectively. Three different possibilities for coefficient A 
and ^ is investigated. The figure demonstrates that the phase diagram 
does not change for local hybridization of impurity with an atom for 
sublattice A or B, or a combination of them. 


the relative weights of the A and ^ changes. In Fig. we 
show comparison with different hybridization patterns. We 
consider three cases: (a) A = 0,^ = 1 (b) A = 1,^ = 0 
and (c) A = 1,^ = 1. The cases (a) and (b) correspond to 
hybridization of the impurity with a molecular orbital of ei¬ 
ther c or d character, while the case (c) above corresponds to 
hybridization with an atomic orbital in sub-lattice A. To have 
a meaningful comparison between the above three cases, we 
should perform an appropriate scaling: From Fermi’s golden 
rule, the broadening of the spin-split impurity levels is pro¬ 
portional to V^. This broadening in (a) and (b) cases is half 
of the (c) case. So we should scale X axis for later case. In 
doing so, the local moment boundaries for all the above three 
cases collapse on the same curve as depicted in Fig.[^ Param¬ 
eter values are indicated in the figure caption. This indicates 
that as long as hybridization remains local, there is no con¬ 
ceptual difference between coupling the impurity to an atom 
from sublattice A, or B or a combination thereof. 


IV. APPLICATION TO MATERIALS AND 
DEFORMATIONS OF THE DIRAC HAMILTONIAN 

As pointed out in the introduction, the derivation by Wolff 
of the anisotropic Dirac Hamiltonian for 3D Dirac solids is 
quite general and applies to a broad range of materials the 
difference of which is reflected in model parameters. The ini¬ 
tial motivation of Wolff was to construct an effective single¬ 
particle theory for the low-lying electronic states of bismuth. 
In this section let us discuss how do the material specific con¬ 
siderations affect the results. We obtain phase diagrams for 
various situations corresponding to variations in different pa¬ 
rameters of the model, and adopting numbers related to bis¬ 
muth. 
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FIG. 4. (Color online) Phase diagram for the local moment formation 
in SIAM for bismuth. The area enclosed by the curves and the Y -axis 
is magnetic region. In this plot the variables X and Y are defined 
by X = \V\^D/{2'kv^) ,Y = (fi — £o)/U. As discussed in the 
text, the top panel shows the hybridization of impurity level with 
valance band the the bottom panel represents the hybridization with 
conduction band. 


corresponding to different values of chemical potential to co¬ 
incide after a vertical shift that compensates the difference in 
the chemical potentials. The fact that after such a shift the 
upper part of the phase diagram coincide indicates that this 
region is basically formed by continuum of states in the va¬ 
lence band. On the other hand, the lower lobe-shaped region 
for three different chemical potentials do not coincide after a 
simple shift and hence they are contributed by the conduction 
band states. This piece of phase diagram is qualitatively close 
to the magnetic region of one-band hosts and 2D graphene. 
Therefore the additional elbow-shaped part of the phase dia¬ 
gram can be considered a feature of two-band 3D Dirac sys¬ 
tems. Having separate contributions to the local moment for¬ 
mation from the two bands of the host material is reminiscent 
of the qualitatively different diamagnetic behavior of bismuth 
compared with normal (one-band, non-Dirac) metals which 
can be understood based on a two-band picture and the 3D 
Dirac HamiltoniarP. The presence of the other band leads to 
two different regions (i) and (ii) in the local moment phase di¬ 
agram of bismuth. This feature is qualitatively different from 
that of normal metals. The lobe-shaped feature of the region 
(ii) is qualitatively similar to the magnetic moment region of 
2D Dirac fermions in graphene, but the elbow-shaped remains 
a feature peculiar to 3D Dirac fermions. 


A. Timing impurity level 

The first thing in a host of bismuth with gap parameter 
A = 7.5 meV can be changed is the type of impurity. This 
obviously changes the value eg. So let us consider the tun¬ 
ing of the parameter Eq of the SIAM. Relative permeability 
of bismuth and its compounds is in the range of 10-40 which 
results in decrease of Coulomb repulsion up to two orders of 
magnitudd^. Following the work of Haldane and Anderson 
on impurity in semiconductor materialP^, we expect that re¬ 
duction of charge accumulation on impurity site to give rise 
to decrease of Coulomb repulsion up to two orders of mag¬ 
nitude. Therefore, we investigate the problem for very small 
values of Coulomb interaction. We focus on the case where 
U = 154 meV { U/A = 20) and present phase diagram for 
three values of chemical potential, p/A = 20,30,40 in Fig.|^ 
This phase diagram has been constructed by varying X and Eq 
( the variation of Y through variation of Eq is emphasized by 
a subscript: Y^^) for three representative values of the chemi¬ 
cal potential p. As can be seen, upon tuning Eq, the magnetic 
region splits into two regions (i) an elbow shaped region for 
larger Y values depicted in the upper part of the figure, and 
(ii) a lobe-shaped region for smaller values of Y. 

We have deliberately separated the regions (i) and (ii) 
above, and have applied a shift S^/U of the whole curves 
along the vertical axis to reveal the different behavior of the 
two regions upon such a vertical shift. Since we have se¬ 
lected the above three values of the chemical potential to be in 
the conduction band, the contribution from integrations over 
the valence band will be identical for three chemical poten¬ 
tials. Therefore we expect different phase boundary curves 


B. Tuning U 

In the single impurity Anderson model for normal metals, 
the combination Y = {^ — eo)/U naturally appears in the 
mean field theory. Therefore as long as the variable Y varies, 
it does not matter which of the three parameters p, Eq, C/ is re- 



X = D/n^ 


FIG. 5. (Color online) Local magnetic moment phase diagram ob¬ 
tained by tuning the Hubbard U. In this diagram the chemical po¬ 
tential and impurity level are set to be 35.4 meV and 20 meV, re¬ 
spectively. By tuning the Coulomb repulsion U in this case where 
fj, — So >0, the magnetic region is qualitatively similar to Fig. [T] 
obtained by tuning fi. 
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derstandable from the atomic limit. If we increase spin-orbit 
coupling further, as already pointed out, we have a region in 
which a small amount of Coulomb repulsion (as small as 5 
meV) leads to formation of local moment. This is probably 
one of the interesting aspects of the local moment formation 
in a three dimensional Dirac material with strong spin-orbit 
interaction that has no counterpart in normal metallic hosts 
without spin-orbit interactions. This means that in a 3D Dirac 
solid even impurity orbitals such as s or p orbitals possess¬ 
ing smaller values of Hubbard U as compared to d or / elec¬ 
tron adatoms get a chance of magnetization! Finally, let us 
focus on the white region adjacent and parallel to Y axis in 
Fig. i that corresponds to the infinitely large spin-orbit cou¬ 
pling, u —?► oo. As the phase diagram shows in this situation, 
irrespective of the value of U, no local moments are formed. 


FIG. 6. (Color online) Local magnetic moment phase diagram ob¬ 
tained by tuning the Hubbard U. In this diagram the impurity level 
lies above the chemical potential p — eo <0. This diagram suggests 
that in this case for Coulomb repulsion as weak as 5 meV, with the aid 
of strong spin orbit coupling, we may find nonzero magnetization. 


sponsible for variation of parameter Y. But since in the case 
of three dimensional Dirac materials the parameter Y does not 
emerge naturally, when constructing the traditional phase dia¬ 
grams in the XY plane the quantity that gives rise to variation 
in Y becomes important. For 3D Dirac solids one natural pa¬ 
rameter is but the other parameters can in principle be 

varied independently leading to a multi-dimensional phase di¬ 
agram. Therefore keeping some fixed, and varying others cor¬ 
responds to viewing a projection of multi-dimensional phase. 
This can be viewed as an advantage as it may reveal new fea¬ 
tures as we will show in this section. Let us see what hap¬ 
pens when we construct the phase diagram in the XY plane 
by tuning the Hubbard parameter U. We consider two cases, 
/i — £0 > 0 and /r — £o < 0. The first case is shown in Fig|^ 
Its general features are similar to Fig. [^obtained by tuning the 
chemical potential p. 

However when we repeat the same analysis for p — eg < 0 
(and U > 0 which means negative values of Y) the phase 
diagram will be completely different and brings about a very 
peculiar physics of local moment formation in 3D Dirac ma¬ 
terials. As can be seen in Fig for a very large values of 
spin-orbit coupling (reflected in v) even with small values 
of U quite strong local magnetic moments can form. This 
corresponds to the red vertical part of the phase diagram in 
Fig. 0 On the contrary looking at the right end of the mag¬ 
netic region in this figure at X ~ 2 axis which corresponds to 
small spin-orbit regime indicates that in this region there are 
no local moment for very small spin-orbit couplings even with 
large Coulomb repulsions. By increasing spin-orbit coupling 
(reducing X), we come to a region in which local moment 
formation is confined between an upper and lower boundary 
values for the Hubbard repulsion. The existence of a mini¬ 
mum repulsion to have unequal populations and is un- 


V. DISCUSSION AND SUMMARY 

Let us first summarize our findings so far; We have in¬ 
vestigated the formation of localized magnetic states in three 
dimensional Dirac solids, and have found that the spin-orbit 
coupling significantly helps with the formation of local mag¬ 
netic moments. Our investigation shows that the effect of 
spin-orbit coupling is to enhance the local moments once they 
are formed. It also allows for formation of local moments 
with very small values of Hubbard U for strong spin-orbit 
couplings. This means that even s or p orbital adatoms whose 
Hubbard U is usually smaller than d or / orbital atoms may 
find a chance of getting magnetized in a host of 3D Dirac 
solid. This chance is not present for them in normal metals. 
Construction of phase diagram by tuning the impurity level 
£o gave rise to two disjoint pieces of magnetic regions. The 
elbow-shaped region having no counterparts in normal met¬ 
als nor in 2D Dirac system (graphene) is due to presence of 
a second band, while the lobe shaped part of the phase dia¬ 
gram comes from the band crossing the Fermi level. In the 
SIAM for a host of three dimensional Dirac material, unlink 
the normal metallic hosts, the parameter Y = {p,—eo)/U does 
not naturally emerge. Therefore in principle the phase dia¬ 
gram should be constructed in a multi-dimensional parameter 
space. Insisting to represent the phase diagrams in terms of 
traditional XY parameters brings in interesting aspects of the 
localized magnetic states in three dimensional Dirac solids. 
Such a larger phase space may provide opportunities for new 
applications and directions in the magnetic properties of of 
Dirac solids. 

Now let us speculate on some deformations of the isotropic 
Wolff Hamiltonian First thing to discuss is the role of 
anisotropy which becomes relevant when it comes to real ma¬ 
terials: In presence of anisotropy, the velocity will be different 
for three different directions, and hence the dispersion relation 
of the host Dirac material becomes, 

Ek = ±^vlkl -F vlkl + vlkl + A2. (17) 

If we use this dispersion in Eq. ( [T3] l when it comes to integra¬ 
tion over k, we can rescale variables as Vxkx —> vkx- This 
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leads to a Jacobian of the form /(vxVyVz) multiplying the 
same integral as the one in the isotropic case. Therefore the 
self-energy rescaled by Jacobian which can be taken in to ac¬ 
count by appropriately redefining v. Therefore, the role of 
anisotropy is just a matter of scaling X axis and does not af¬ 
fect the qualitative physics discussed in this paper. 

There is another conceptually important deformation of the 
simple Dirac Hamiltonian. The Dirac Hamiltonian Q can 
also be generalized by adding a quadratic term 

3 

Hg = vY^ + {A- BP)(3. (18) 

where and 13 = ^^ are usually defined in Dirac 

equation in terms of Dirac matrices 7 ^. This generalization 
allows for two possibilities with respect to topology of the re¬ 
sulting host material: When BA < 0 it is topologically trivial, 
while if BA > 0 it has a non-trivial topology^ with associ¬ 
ated boundary states. Having gained some insight into the 


important role of spin-orbit coupling in the local moment for¬ 
mation in 3D Dirac solids, we can briefly address the role of 
quadratic B term in the limit of small B. In this Taylor expan¬ 
sion of the resulting dispersion relation leads to a straightfor¬ 
ward renormalization of velocity i.e. ^ — BA. There¬ 

fore, in topologically non-trivial (trivial) case where BA > 0 
{BA < 0), the quadratic term leads to a decrease (increase) 
in the effective spin-orbit interaction. Therefore the topolog¬ 
ical twist of spins in the momentum space corresponding to 
non-trivial topology are expected to weaken the aspects of lo¬ 
cal moment physics of 3D Dirac solids discussed in this pa¬ 
per, while in the topologically trivial case, at least within the 
present perturbative scheme limited to very small B, the spin- 
orbit driven aspects of local moment physics are expected to 
get enhanced upon addition of the quadratic term B. 
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